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We perform a low temperature Raman scattering study of phononic and collective spin excita- 
tions in the orthogonal dimers compound SrCu2(B03)2 , focussing on the symmetry and the effects 
of external fields on the magnetic modes. The zero field symmetry and the behavior in magnetic 
fields of the elementary and bound magnetic triplet states are experimentally determined. We find 
that a minimal 4-spin cluster forming the unit cell is able to describe the symmetry as well as 
the anisotropic dispersions in external fields of the spin gap multiplet branches around 24 cm -1 . 
We identify two Raman coupling mechanisms responsible for the distinct resonance behavior of 
these magnetic modes and we show that one of these can be ascribed to an effective intra-dimer 
Dzyaloshinskii-Moriya spin interaction. Our data also suggest a possible explanation for the ex- 
istence of a strongly bound two-triplet state in the singlet sector which has an energy below the 
spin gap. The low temperature phononic spectra suggest strong spin-phonon coupling and show 
intriguing quasi-degeneracy of modes in the context of the present crystal structure determination. 

PACS numbers: 75.10.Jm, 75.50.Ee, 78.30.-j, 71.70.Ej 



I. INTRODUCTION 

Several properties make SrCu2(B03)2 a unique system 
among the known quantum magnets. This compound 
is a realization of a two-dimensional (2D) spin system 
with a disordered ground state even at very low temper- 
atures and a spin gap A « 24 cm -1 (3 meV) from the 
singlet (S — 0) ground state to the lowest excited triplet 
(S = 1) stated. The strengths of the relevant magnetic 
interactions place this compound close to a quantum crit- 
ical point (QCP). Moreover, data in high magnetic fields 
show plateaus at commensurate (1/8, 1/4 and 1 /3) values 
of the saturation magnetization^* 3 ^. The plateau states 
can be thought of as crystalline arrangements of mag- 
netic moments separating regions of continuous rise in 
magnetization, the latter allowing for an interpretation 
in terms of Bose-Einstein condensation of triplet excita- 
tions^. It has also been suggested^ that doping in this 
system (regarded as a Mott-Hubbard insulator) may lead 
to a superconducting phase mediated by antiferromag- 
netic (AF) fluctuations, a mechanism similar in spirit to 
one of the scenarios proposed for the high T c cuprates 7 . 

The S = 1/2 Cu spins are arranged in weakly cou- 
pled 2D layers defining the (ab) plane. In each of these 
sheets they form orthogonal spin-dimer lattices^, see 
Fig. la. At T c = 395 K the system undergoes a second 
order phase transition from the space group 1 4/ 'mem to 
/42m on cooling down from the high temperature side. 
In the 1 4/ mem phase the planes containing the Cu atoms 
are flat and they form mirror symmetry elements. The 
transition at 395 K can be understood as the buckling 
of the Cu planes which lose their mirror symmetry prop- 
erty. Inversion symmetry is lost as well below T c but the 
number of atoms in the unit cell remains unchanged due 
to the orthogonality of the spin-dimer network. Because 



of the same orthogonal arrangement, there are two spin 
dimers in the unit cell. As a result, the spin gap exci- 
tation, defined as the transition from the S — ground 
state to that (usually the lowest in energy) S — 1 level 
whose wavefunction contains a single spin dimer excited 
to the triplet stataii2*ii, has a 'fine structure' made out 
of six levels, three from each dimer in the unit cell. 

The 2D spin-dimer lattice can be described well by 
taking into account the nearest and next nearest neighbor 




FIG. 1: (Color online) (a) The structure of the spin-dimer 
compound SrCu2(BOs)2 . The unit cell is shown by the short 
dashed square. The circles are the 5 = 1/2 Cu spins, the 
solid and long dashed lines represent the intra-dimer (Ji) and 
inter-dimer (J2) AF superexchange interactions. The arrows 
correspond to the proposed intra-dimer antisymmetric spin 
interaction d a ;,. Below 395 K the orthogonal dimers are no 
longer coplanar. (b) The 4-spin cluster corresponding to the 
unit cell in (a). As expected, the symmetry of this object 
is D 2d, which is the point group associated with the crys- 
tal space group of SrCu2(BOa)2 below 395 K. Note that the 
inter-dimer superexchange J2 is effectively doubled if periodic 
boundary conditions are used. 
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AF super-exchange interactions J\ and Ji respectively, 
see Fig. la. In this approximation and for a ratio x — 
J2/J1 lower than about 0.7, the direct product of singlet 
dimers on each rung, which is always an exact eigenstate 
of the Hamiltonian, is also the ground state of the system 
and the spin gap is finite although it gets renormalizcd 
down with increasing x due to many body effectsAfl^. 
For large x the ground state changes and the Hamiltonian 
has long range AF order, other possible intervening states 
separated by QCP's being proposed to exist around 0.7, 
see Refi Based on fits to magnetization data, predictions 
of symmetries, relative energies and dispersions of single 
and composite triplet excitations, theoretical estimates 
for x range from 0.603 Refsi2i±i to about 0.68 Refsi 2 -^ 

An analysis of the vibrational modes is of interest in 
SrCu2(BC>3) 2 because spin-lattice interaction has been 
suggested to be relevant to the magnetic dynamics at 
low temperatures and high magnetic fields. In particular, 
spin-phonon interaction has been invoked in order to ex- 
plain the selection rules of the magnetic transitions seen 
in infra-red (IR) absorption^. The coupling between the 
lattice and magnetic degrees of freedom was also taken 
into account in order to describe the spin density profile 
at high fields in the magnetization plateau states^. So 
far the study of phononic excitations has been related 
mostly to the crystallographic changes at 395 K. It has 
been established in RefsAi^ that the structural transi- 
tion displays soft mode behavior (see RefsAi^ also for a 
phononic symmetry analysis). The phononic mode which 
condenses at 395 K belongs to the Bi„ representation of 
the high temperature symmetry group and, in terms of 
Cu atoms, it involves essentially an alternate displace- 
ment along the c-axis of the nearest neighbor dimers. In- 
terestingly, several phonons appearing below T c were ob- 
served in (ca) polarization as very close in energy shoul- 
ders of some of the modes found above the transition. 
Due to the absence of inversion symmetry below 395 K 
this experimental fact suggested the existence of almost 
degenerate even and odd (transforming like the E g and 
E„ representations respectively) excitations in the high 
temperature phaseAa. 

In terms of magnetic properties, elementary and two- 
triplet excitations have been studied by inelastic neutron 
scattering (INS)i£»A£, electron spin resonance (ESR)i2iiS 
and IR absorption spectroscopy 15 20 . The low temper- 
ature Raman scattering results shown in Refi^i are fo- 
cussed on composite magnetic excitations in the 5 = 
channel and are used to extract quantitatively the ratio 
x = J2 1 J\ ■ In our Raman study, the emphasis is however 
on the elementary and composite collective excitations 
in the 5 = 1 sector, the main results being related to 
the experimental determination of their symmetries and 
anisotropic behavior in magnetic fields. We are also able 
to relate some of these new findings to results of a 4-spin 
cluster analysis as well as to propose an additional effec- 
tive spin interaction induced by spin-orbit coupling. It 
has been established that due to the frustrated nature 
of the magnetic interactions, the one triplet excitations 



are local, weakly dispersive in the reciprocal space, while 
two-particle states are more mobile^ and have contribu- 
tions from the whole Brillouin zone 1011 . 

In spite of a lot of experimental and theoretical ef- 
fort for understanding the magnetic properties there are 
several open questions. An exact determination of the 
ratio of the exchange interactions x, which is important 
due to the proximity to the QCP, is still missing. One 
interesting aspect in this regard is the observation of a 
magnetic state at 21.5 cm -1 , which is below the energy of 
the spin gap multiplet A. This could be seen in ESR 19 , 
INSii and IR™ data in energy level anti-crossings in the 
downward dispersions with magnetic field of some of the 
spin gap branches. The existence of such a low energy 
mode brings into question the exact quantitative estima- 
tion of the AF exchange parameters. As we will show, a 
direct comparison between our experimental findings and 
theoretical predictions regarding the symmetry of these 
modes is illuminating in this respect. 

A different set of questions is related to the way the 
external radiation field couples to the magnetic excita- 
tions. While the photon induced spin-exchange process 
insures the Raman coupling to 5 = two-triplet exci- 
tationsSii 2 ^, the transitions to 5 = 1 states require the 
presence of spin-orbit coupling 22 . Although an effective 
antisymmetric Dzyaloshinskii-Moriya (DM) term origi- 
nating in the spin-orbit interaction (the contribu- 
tion explicitly written in Eq. has been proposed to 
explain the 'fine structure' of the six levels forming the 
gap multiplet around 24 cm -1 seen in a neutron scat- 
tering study, Ref^i, this term does not mix the singlet 
ground state to excited 5 = 1 states. As a result, the na- 
ture of the coupling leading to the observed ESR and IR 
data is still to be understood. Possible candidates were 
suggested in the discussion of ESR datai^ and dynami- 
cal spin-phonon induced effective DM interactions have 
been invoked in order to explain the optical absorption 

In this article we study phononic modes and collec- 
tive magnetic excitations in SrCu2 (603)2 with the ex- 
perimental emphasis placed on the understanding of the 
magnetic Raman scattering in the 5=1 channel. In the 
low temperature in-plane polarized phononic spectra we 
find (in the 50 - 350 cm -1 energy range) several pairs of 
modes with similar energies. This quasi-degeneracy is in- 
triguing because group theory analysis suggests that they 
are related to different atomic vibrational patterns. Re- 
garding the magnetic Raman scattering data, our novel 
findings are the following. The symmetry analysis and 
exact diagonalization of the 4-spin system Hamiltonian 
(see Fig. lb) given by: 

h = h + h? m + h3, m (1) 

which includes the main, unperturbed, Heisenberg term, 
H , the antisymmetric inter-dimer DM term proposed in 
Refiii, H® M , and an additional intra-dimer DM term, 
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H° M , reading: 

Ho = J\ Si • Sj + J2 E] Si ■ Sj, (2) 

#-? M = E d c (!j) (S t xS 3 ), (3) 

(ij) NNN 

HT= E ^(SiXSi) (4) 

(i,3) WW 

explains the experimentally determined symmetries of 
the zero field Brillouin zone center spin gap branches 
around 24 cm -1 (confirming the local nature of the ele- 
mentary one-triplet modes) but fails to account for the 
two-triplet excitations. In the equations above i and j are 
nearest neighbor (NN) or next nearest neighbor (NNN) 
Cu sites and d c and d a b are inter and intra-dimer anti- 
symmetric spin interaction vectors. Remaining confined 
to the 4-spin cluster we find that by introducing the effec- 
tive intra-dimer DM interaction H® b M we are also able 
to reproduce the observed selection rules and intensity 
variations of the spin gap branches in external magnetic 
fields. These selection rules also require that the energy 
of the S = two-triplet bound state formed from spins 
confined within a unit cell is below A (in the 4-spin clus- 
ter this is equivalent to x > 0.5 in Fig. 4) suggesting a 
high binding energy for this two particle excitation. In 
the last section, we also show results of a resonance study 
which allows us to identify the action of two different light 
coupling mechanisms to magnetic excitations. 

II. EXPERIMENTAL 

Raman scattering from the ab surface of a single crystal 
of SrCu2 (1303)2 , grown as described in RefiS^, was per- 
formed using an incident power density less than 1 mW 
focussed to a 100 /j,m diameter spot. The crystallographic 
axes orientation was determined by X-ray diffraction. 
The data in magnetic fields, taken at a sample tempera- 
ture of about 3 K, was acquired with a continuous flow 
cryostat inserted in the horizontal bore of a supercon- 
ducting magnet. We used the ul — 1-92 and 2.6 eV 
excitation energies of a Kr + laser and a triple-grating 
spectrometer for the analysis of the scattered light. The 
data were corrected for the spectral response of the spec- 
trometer and CCD detector. 

Polarized Raman scattering can probe Brillouin zone 
center excitations that belong to different symmetry rep- 
resentations within the space group of the crystal struc- 
ture. We denote by (ej„e ottt ) a configuration in which 
the incoming/outgoing photons are polarized along the 
e»n/ e ut directions (see Fig. 1 for axes notations). The 
(RR) and (RL) notations refer to circular polarizations, 
e in = (a - ib)/V2, with e out = e in for (RR) and 
Gout — e i n for the (RL) geometry. The point group as- 
sociated to the high temperature phase space group of 
the SrCu2(B03)2 crystal, IA/mcm, is D^. The point 



group associated to the /42m space group corresponding 
to the lower temperature phase is Z?2d- I n ^2di the (RR), 
(RL), (aa), (ab), (a 1 a 1 ) and (a'b 1 ) polarizations probe ex- 
citations which belong to the Ai + A2, Bi + B2, Ai + 
Bi, A 2 + B 2 , Ai + B 2 and A 2 + Bi irreducible repre- 
sentations of D2d- 

III. PHONONIC SYMMETRIES 

In Fig. 2 we show six low temperature Raman spectra 
excited with the u>l = 1-92 eV laser energy. The arrows 
point to the observed modes above 60 cm -1 and also 
shown arc their energies and symmetries. Magnetic fields 
do not affect their energies which indicates the phononic 
nature of these excitations. The modes below 60 cm -1 
are not indexed since they will be discussed in the next 
sections. The excitation at 59 cm" 1 corresponds to the 
fully symmetric soft mode of the structural transition at 
395 K^iil. The corresponding two-phonon excitation in 
the Ai channel is seen at 121.8 cm -1 and very close to it 
another sharp mode with Bi symmetry. Similarly to the 
data in (ca) polarization in Refr 4 , where below 395 K 
several new modes appear as shoulders of phonons ex- 
isting above the transition, we observe several pairs of 
modes having close energies. For example doublet struc- 
tures are observed around 155 cm -1 , where we see a pair 
of B2 and A2 excitations, around 285 cm -1 , where we 
observe two modes having Bi and Ai symmetries, and 
also around 320 cm -1 where we see a pair made of B2 
and Bi symmetric excitations. 

One way to explain this behavior, suggested in Refii^, 
is to assume that in the high temperature phase there are 
phonons which are odd (ungerade: V) and even (gerade: 
'g') with respect to inversion but very close in energy and 
to try to identify them by looking at similar atomic vi- 
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FIG. 2: (Color online) Raman spectra at T = 3 K in six po- 
larization configurations and zero applied field. Next to each 
vertical arrow we show the energies (in cm -1 ) and symmetries 
of the modes above ~ 60 cm -1 . The spectra are vertically off- 
set. The excitation energy used is u>l = 1-92 eV. 
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brations corresponding to 'u' and 'g' representations re- 
spectively. Following this idea in more detail and taking 
into account that the set of irreducible representations 
{Ai g ,A2g,Bi ff ,B 2g ,Ai M ,A2„,Bi M ,B2n} of the D ih point 
group becomes the set of {Ai,A 2 ,Bi,B 2 ,Bi,B2,Ai,A2} 
representations (in this order) in Z?2d, we would have 
for example that the (Ai,Bi) pair around 284 cm -1 cor- 
responds either to a (Ai ff ,Ai u ) group or to a (Bi„,Bi g ) 
group in the high temperature phase. This is because 
we chose gerade-ungerade pairs and, as shown above, the 
Ai symmetric mode can originate either from a A lg or 
a Bi M symmetric phonon while the Bi mode could be 
either a Bi g or a Ai u phonon in the high temperature 
phase. Similar reasoning would suggest that the (Bi,B2) 
group around 320 cm^ 1 originates either from a pair of 
(Bi fl ,A2u) or (B2 ff ,Ai u ) above 395 K and also that the ori- 
gin of the (B2,A2) group around 155 cm -1 can be pairs 
of (B2 5 ,B2 M ) or (A2 g ,Bi„) modes in the 1 4/ mem phase 
which has as the associated point group. 

We performed a symmetry analysis of the k = atomic 
vibrations in the high temperature phase and our con- 
clusion is that the approach suggested in Ref. 14 does not 
provide an a priori reason for the quasi-degeneracy. This 
conclusion, as explained in the following, is based just 
on a simple inspection of the character table of the 
point groupSa. One can easily note that in D±h the even 
modes are symmetric with respect to the mirror symme- 
try in the Cu(B03) planes while the odd vibrations are 
antisymmetric with respect to this symmetry operation. 
This means that the 'u' phonons in the high tempera- 
ture phase correspond to vibrations of the atoms along 
the c-axis while the 'g' modes consist of in-plane move- 
ments. Due to this difference in the oscillation patterns, 
one parallel and one perpendicular to the CUBO3 planes, 
it is hard to explain the closeness of phononic energies at 
this qualitative level. This is why we find intriguing the 
observed quasi-degeneracy of the phononic modes in the 
context of the present crystal structure determination. 

We remark that on the other hand one could find vi- 
brations which involve similar oscillations at the 'molec- 
ular' level (for instance groups of atoms forming the Cu- 
O spin dimer structure or groups of 02 atoms bridging 
nearest neighbor spin dimers) and which belong to differ- 
ent group representations because of the different 'inter- 
molecular' phase pattern. We suggest that, remaining 
within the conclusions of X-ray studies^ which so far 
have not found evidences for additional crystallographic 
changes at low temperatures (that in turn may produce 
phonon splittings), good candidates for understanding 
this behavior are provided by the inter-dimer BO3 molec- 
ular complexes whose rotations as a whole around the a 1 
and c-axes or whose in and out of the plane translations 
may turn out to be similar in energies. The data suggest 
that theoretical work in this respect could be interesting. 

We note that the appearance of the weak 155.9 cm -1 
mode in the (RR), (ab) and ia'b') polarizations singles 
this mode out from the other excitations because it be- 
longs to the A2 symmetry channel. The Raman cou- 
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FIG. 3: (Color online) Zero field Raman data in taken with 
u>l = 1-92 eV excitation energy at T = 3 K in six polarization 
configurations. The legend shows the tetragonal symmetries 
probed in each scattering geometry. See the text for nota- 
tions. 

pling to this excitation is unusual compared to the other 
modes in the sense that it cannot take place via two 
electric dipole transitions. We find several magnetic res- 
onances (see Fig. 3) with A2 symmetry, but the absence 
of magnetic field effects suggests that this mode has a 
preponderant phononic character. A finite intensity of 
such excitation provides direct evidence for spin-phonon 
coupling in SrCu2(B03)2 • 



IV. MAGNETIC SYMMETRIES 

The irreducible representations of the I?2d point group 
probed by six scattering geometries, discussed also in the 
experimental section, are shown in the legend of Fig. 3. 
In this figure, the six low temperature Raman spectra 
from Fig. 2 are shown for the frequency region below 
60 cm -1 . Three strong features around 23, 29 and 38 
cm" 1 , denoted by T^,, Si and T2, are observed and they 
transform like the A2, Bi and A2 representations respec- 
tively. Besides these three modes, we observe several 
other weaker excitations. Among them we see a set of 
three A2 symmetric modes denoted by T3, T4 and T5. 
We also note the presence of the excitations denoted by 
Ti e and Tiy which have Bi and B2 symmetries, giving 
rise to a small feature seen around 25.6 cm -1 in all po- 
larizations except for the (RR) configuration. Table [I] 
contains a summary of the observed excitations below 
the 60 cm -1 frequency range. 

The energies of these Raman active excitations are in 
agreement with those where INSi&, ESRi&ia and IRiS^ 
data observed magnetic modes. The data in magnetic 
fields (see the discussion in the following section) con- 
firm the magnetic nature and the predominant S = 1 
character of the 'T' modes. Therefore T^,, Ti e and Tif 
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FIG. 4: (Color online) (a) The eigenvalues of the Ho term 
of the spin Hamiltonian of Eq. corresponding to the 4-spin 
cluster shown in Fig. lb. The spin gap excitation at A = Ji 
is six times degenerate in this approximation, (b) Symmetry 
analysis of the 4-spin cluster in (a) in the Dm group. The 
energies on the left are calculated for a ratio x — J 2 / Ji equal 
to 0.3. On the right we show the energy splittings and the 
absolute and relative to the ground state (in parenthesis) sym- 
metries of the 16 magnetic modes when a finite inter-dimer 
DM interaction d c (see Ref^-) is present. The effect of a fi- 
nite intra-dimer DM perturbation, d a j,, is to further split the 
modes which belong to the one-dimensional representations. 



excitations modes seem to belong to the spin gap multi- 
plet while T2, T3, T4 and T5 would correspond to multi- 
particle triplet channels. 

If the picture of real space localized elementary triplets 
is true, then one expects that the analysis of the 4-spin 
cluster in Fig. lb forming the unit cell is able to pre- 
dict correctly at least some of the experimentally ob- 
served symmetries of these excitations. Besides symme- 
try analysis, numerical diagonalization of the Hamilto- 
nian of Eq. allows one in principle to also identify the 
energy as well as the predominant spin character of each 
eigenstate of this cluster. Indeed, as can be seen by com- 
paring Figs. 3 and 4b and looking at the rightmost col- 
umn of the Table [J the symmetries of the observed zero 
field one-triplet excitations correspond to the results of 
group theory analysis. The doubly degenerate E modes, 



TABLE I: Collective spin excitations in zero field: notation, 
the predominant spin character (Stot), the z projection of the 
spin (S z ), the energies and transition symmetries as observed 
experimentally and predicted from the 4-spin cluster in Fig. 1 
corresponding to k — excitations. T\ a ... Tif represent 
elementary triplet excitations. Generally, the modes whose 
energies change in external magnetic fields are indexed by T. 
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Experiment 


Group Theory 




1 (±l) 


22.8 




Ai 


Tib 


1 (±l) 


22.8 


A 2 


A 2 


Tic 


1 (0) 


24.2 
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T w 


1 (0) 


24.2 




E 


Tie 


1 (±1) 


25.6 


Bi 


Bi 


Ti/ 


1 (±1) 


25.6 


B 2 


B 2 


Si 


0(0) 


28.9 


Bi 




T 2 


1 


37.5 


A 2 




T 3 


1 


40.8 


A 2 




S 2 





43.0 


Bi 




T 4 


1 


44.5 


A 2 




S3 





45.2 


Bi 




T 5 


1 


50.9 


A 2 





Tic and T\d, are not observed in zero field when the light 
propagates parallel to the c-axis because they are acces- 
sible only in (ca) or (cb) polarizations. The observation 
of the Ti e (Ti/) modes with Bi (B2) symmetries at an 
energy of 2.8 cm -1 above the A2 symmetric mode Tib 
allows the determination of the magnitude of the inter- 
dimer interaction d c (see Fig. 4b) and also of its absolute 
signSi. The fully symmetric Ti a mode which, within the 
spin model including only the Ho and H® M terms given 
in Eqs. (0 and ©, should be degenerate with the strong 
Ti;, (A-2 symmetric) excitation at 22.8 cm -1 is also not 
observed which is most probably due to a much weaker 
coupling to light in this symmetry channel. We will dis- 
cuss the coupling mechanisms in the last section of our 
paper. 

We turn now to the discussion of two-triplet states. Be- 
sides the singlet ground state and six one-triplet states, as 
can be seen in Fig. 4, the 4-spin cluster analysis predicts 
the following: one S — two-triplet bound state in the 
Bi channel, three branches with Ai and doubly degener- 
ate E symmetries which belong to another bound S = 1 
excitation and five branches of a quintuplet (S — 2) state 
having Ai, A2, Bi and E symmetries with respect to the 
ground state. Their energies are plotted in Fig. 4a as a 
function of x = Jij J\. We note that due to symmetry 
reasons none of the observed A2 symmetric modes from 
T2 to T5, having energies higher than 30 cm" 1 , qualify 
for an interpretation as triplet bound states generated 
within the 4-spin cluster. This is consistent with the fact 
that larger cluster sizes are necessary in order to capture 
the more delocalized nature of these excitations. 

The fact that the existence of the strong A2 symmet- 
ric bound triplet state at an energy 1.55 -A = 37.5 cm" 1 
has not been predicted by high order perturbative anal- 
ysis^ (here we refer especially to the symmetry of this 
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excitation, not its energy) suggests that other spin inter- 
actions have to be taken into account in order to explain 
the excitation spectrum. Apparently the symmetry con- 
siderations would allow the 28.9 cm -1 feature denoted 
by Si in Fig. 3 to be interpreted as the singlet bound 
state of two triplets within a unit cell. As we show in 
the following section, the 28.9 cm -1 mode does not shift 
in external fields, which is compatible with a collective 
singlet excitation as discussed in2i, but suggests that its 
internal structure is not the one derived from the 4-spin 
cluster. 



V. MAGNETIC FIELD EFFECTS AT T rs 3 K 

In Fig. 5, using the same mode notations, we show the 
influence of an external magnetic field applied parallel 
and perpendicular to the c-axis on the low temperature 
Raman spectra from Fig. 3. Here we summarize the rel- 
evant aspects, noting first that energy shifts induced by 
magnetic fields in this figure were observed only for the 
modes indexed by T in the Table [I] In Fig. 5a we ob- 
serve the splitting of the Ti a and Tib modes in magnetic 
fields B || c, the B = 1 T spectrum showing that the A 2 
mode (Tib) present in zero field disperses upwards with 
increasing magnetic field. Dashed lines mark the disper- 
sion of the much weaker modes T3 , T4 and T5 . In Fig. 5b 
one of the E modes at 24.2 cm -1 becomes Raman active 
due to symmetry lowering for B _L c configuration and 
we observe three dispersing branches of the gap multi- 
plet. Fig. 5c shows that the Bi symmetric excitation at 
28.9 cm -1 does not change its energy with field, only a 
very small negative shift of the order of 0.5 cm -1 from 
to 6 T is seen because of the crossing with the up- 
ward dispersing gap branches seen in (RR) polarization. 
Fig. 5d, which is a zoomed in region of Fig. 5b, shows 
that several modes become Raman active in finite fields 
B _L c around 38 cm -1 where the T2 excitation lies. The 
internal structure of this higher energy multiplet is com- 
posed of modes dispersing up, down and independent of 
magnetic field. We remark on the similarity in selection 
rules and dynamics in magnetic fields between the col- 
lective modes around 38 cm -1 and that of the spin gap 
branches around 24 cm . The emergence in finite fields 
of several strong modes in the spin gap region precludes 
the observation of the dynamics of the weak Ti e and Tiy 
modes from Fig. 3. 

Fig. 6 summarizes the magnetic field dependencies of 
the energies and spectral weights of the most intense Ra- 
man excitations. The symbols in Fig. 6a-d correspond 
to the experimental data and the solid lines are results 
of calculations: the energies in panels a and b are ob- 
tained by exact diagonalization of the Hamiltonian of 
Eq. ^ with the parameters specified in the next para- 
graph; using Fermi's golden rule, the intensities in Fig. 6c 
and Fig. 6d are calculated as the square of the matrix ele- 
ments between the ground and excited states of the effec- 
tive Fleury-Loudon spin interactions describing the cou- 
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FIG. 5: (Color online) Magnetic field dependencies of the 
magnetic excitations at T = 3 K using the ujl = 1.92 eV 
excitation energy in the following geometries: (a) (RR) B \\ c. 
(b) (RR) B _L c. In (c) the (RL) polarized data is shown for 
and 6 T magnetic fields for both B \\ c and B _L c. In (a) 
and (b) the vertical shift is proportional to the magnetic field 
difference with respect to the 9 T spectrum and the dashed 
lines are guides for the eye. Panel (d) is a zoomed in region 
of panel (b). 



pling to the external electromagnetic field 22 . The form of 
these interaction terms, denoted by iJ™* and H™*, are 
discussed explicitly in the next section. 

Here we discuss the parameter set used for data fitting 
in both the B \\ c and B _L c configurations. Taking 
into account that: (i) the 4-spin cluster neglects many- 
body gap renormalization effects leading to a singlet- 
triplet energy independent of x = J2/J1, see Fig. 4a, 
and (ii) the fact that when using periodic boundary con- 
ditions there is an effective doubling of the J2 and inter- 
dimer DM interactions, we chose the following values: 
Ji = A = 24.2 cm -1 which is the value of the spin gap 
(see the Table and Refsi 16 i 17 i 18 i 19 i 20 ); 

x = 0.556 from the ratio of the energies of the sub-gap 
mode at 21.5 cm" 1 , Refs4&*i£, with respect to the gap A 
(see Fig. 4a); an inter-dimer DM term d c = 1.4 cm -1 , 
which produces the splitting of the Ti a b and Tiej 
branches from 24.2 cm -1 to 22.8 and 25.5 cm -1 respec- 
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FIG. 6: (Color online) Energies (panels a and b) and intensi- 
ties (panels c and d) of the spin excitations for B || c (left) and 
B _L c (right) from Fig. 5a-b. Symbols represent experimental 
points, solid lines are the results of 4-spin cluster diagonaliza- 
tion as described in the text, dashed lines are guides for the 
eye. 



tively (our value of d c is consistent with the one pro- 
posed in Ref^i); finally, using the magnetic field value 
where the intensities cross in Fig. 6c, the value of the 
intra-dimer interaction was set to d a b — 2.66 cm -1 . We 
note two aspects regarding the magnitude of the param- 
eters used above. The first is that because many-body 
spin interactions are not captured within this minimal 
4-spin cluster, our chosen value for x = J%l J\ should not 
be taken ad litteram for the real structure. The second 
aspect, discussed in more detail later in this section in 
connection to the dispersion of the magnetic modes in 
the B || c configuration, is related to the fact that we 
were forced to choose x > 0.5, which is equivalent to say- 
ing that the singlet bound state has an energy which is 
lower than the elementary gap excitation at 24.2 cm -1 . 
This fact may be important regarding the experimental 
observation of a mode at 21.5 cm -1 , see Refsii^i^iiS. 

The space group symmetry of the crystal uniquely im- 
poses for any existing static intra-dimer DM interaction 
(the term H® b M in Eq. the configuration depicted in 
Fig. la, i.e. the DM vectors are perpendicular to the o 
axis and dimer bonds. It is the H® b M term in the system 
Hamiltonian which is responsible in our interpretation for 
the mixing of singlet and triplet modes and allows for a fi- 
nite coupling of the latter excitations to the external pho- 
ton field. As we mentioned in the introduction, a finite 
intra-dimer interaction d a b is crucial because the H® M 
term alone does not mix the singlet and triplet states. 
Also shown in Fig. 6 by filled squares and dashed lines 
are the experimental field dependencies of other higher 



energy modes observed in Fig. 4. 

The agreement for this set of parameters is qualita- 
tively good overall and quantitatively better with regard 
to the energy and intensity variations for the B _L c case. 
As discussed before, the term H® b M plays a crucial role in 
obtaining a finite coupling to the excited S — 1 triplets. 
H® b M also produces splittings of the Ti a and Tu, modes, 
of the Ti e and Ti / modes as well as that of the quin- 
tuplet branch (these splittings are not shown explicitly 
in Fig. 4b). For the spin gap branches, the magnitude of 
the splittings is unresolved because it is very small, of the 
order of d a 2 b /A w 0.25 cm -1 . The largest discrepancy be- 
tween the experimental data and the calculation is seen 
in Fig. 6c. One aspect in this regard is that the value 
d a b had to be chosen greater than that of d c . This is un- 
expected because d c is allowed by symmetry both above 
and below the structural phase transition at 395 K while 
the existence of a finite intra-dimer DM interaction is al- 
lowed only below 395 K, when the mirror symmetry of the 
(ab) plane is broken. Additional terms may be responsi- 
ble for this disagreement, possible candidates being the 
in-plane components of the inter-dimer DM interaction, 
which should also be allowed below the structural phase 
transition. 

We now discuss the existence of a magnetic mode be- 
low the spin gap valuei^iiiS. In order to reproduce the 
upward dispersion with fields B \\ c of the Tif, mode 
we had to choose a value for x which is greater than 
0.5, otherwise this excitation would have displayed a 
downward dispersion. ^From Fig. 4a we observe that 
x = J^l J\ > 0.5 implies that ^the position of the bound 
singlet state is below A. We suggest that precisely 
this state may be responsible for the observations of the 
21.5 cm -1 mode in Refsii^*i2*iS. The presence of this 
excitation will also influence specific heat measurements 
and, in conjunction with the finite intra-dimer interac- 
tion d a b, also the low temperature magnetization data 
which is not quantitatively understood yeli 

The existence of this magnetic mode below the gap 
seems quantitatively at odds with theoretical predic- 
|.j on9 io 4 i2 i i3^ Nevertheless, perturbational calculations 
predict the existence of a singlet state at 25 cm -1 , which 
is above but very close to the spin gapiS. A quantita- 
tive reconciliation between theory and the observed se- 
lection rules in magnetic fields could be achieved if the 
coupling ratio x = J2/J1 is slightly increased from the 
value of 0.603, as determined in RefpiS. The experimen- 
tal finding of the set of A2 modes (T2 to T5 in Fig. 3, 
all of them below the two-magnon continuum starting at 
2 A ps 48 cm -1 and whose symmetries are also not pre- 
dicted by theory) shows that although several aspects of 
the magnetic bound states are understood, a complete 
picture of the spin dynamics in the multi-triplet sectors 
is still to be achieved. 
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FIG. 7: (Color online) T = 3 K Raman data in (a'b') po- 
larization for ujl = 1-92 (filled squares) and 2.6 eV (empty 
circles) excitation energies. The data points corresponding to 
the resonantly enhanced mode at 29 cm -1 are multiplied by 
0.15 and the corresponding line represents a Lorentzian fit. 



VI. RESONANCE AND LIGHT COUPLING 
MECHANISMS 

Fig. 7 shows two low temperature Raman spectra taken 
in (a'b') polarization with two incoming laser frequencies, 
u)l = 1.92 and 2.6 eV. The point we want to make in 
this paragraph is that we observe two qualitatively dif- 
ferent behaviors. Firstly we notice that the area under 
the peak corresponding to the Tif, mode in the spectrum 
taken with the excitation energy ll>l = 1-92 is of the same 
order of magnitude as the one in the spectrum taken with 
the 2.6 eV laser excitation frequency. The same obser- 
vation is also true for the mode denoted by T2. ^From 
this perspective, a different behavior is observed for other 
^modes and an example is the group of modes denoted 
by Ti e , iSx, S2 and S3. In this latter case for instance, 
it is clearly seen in our spectra that the mode Si in the 
spectrum taken with the lol = 2.6 eV excitation is sev- 
eral orders of magnitude stronger than when using the 
uil = 1-92 eV laser line. Also, while the modes T ie , S2 
and S3 are barely visible when ujl = 1-92 eV (see also 
Fig. 3), they become quite strong for u>l = 2.6 eV. Fig. 5 
shows in addition that other modes in the 50 to 70 cm -1 
energy range also become visible only in the lol = 2.6 eV 
spectrum. We conclude this paragraph by saying that the 
group {Tib,T 2 } has a different resonance behavior than 
the group {Ti e , Si, 82,83} in the sense that their relative 
intensities are quite different for the two laser excitation 
frequencies used. 

Note that in the paragraph above we commented only 
on the relative intensities of the Raman modes when the 
u>l = 1.92 eV excitation was used compared to the 2.6 eV 
line. Although the absolute intensities of each the Raman 



excitations can be affected by optical corrections of the 
data, this effect cancels out when relative intensity ratios 
arc considered. In other words, it is possible in principle 
that optical corrections can change the ratio of the ab- 
solute values of the Raman intensities of the Tu or T 2 
modes when taken with u>l — 1-92 and 2.6 eV laser en- 
ergies respectively. The same observation applies when 
we relatively compare the intensities of the Ti e , Si, S2 
and S3 for these two excitations. However, the important 
point is that the ratio of intensities should be affected in 
exactly the same way if the same resonance mechanism 
was responsible for both the { Tib, T%) group of excitations 
on one hand and for the group of modes {T\ e) 81,82, S3} 
on the other hand. The conclusion following from the 
above discussion is that the data in Fig. 7 proves the 
existence of two light coupling mechanisms to the ob- 
served low energy excitations. In fact, further support for 
our statement comes from preliminary data which show 
that the 1.92 and 2.6 eV laser excitations energies corre- 
spond to regions in the visible range where distinct fea- 
tures of the reflectivity spectrum are observed^. While 
certainly interesting in this respect, a complete determi- 
nation of the optical parameters of SrCu2(B03)2 in the 
visible range and the identification of the specific high 
energy electronic intermediate states involved in the cou- 
pling mechanisms goes beyond the scope of our article 
and, as explained above, cannot change our conclusions. 

The resonance of the Bi symmetric T e magnetic mode, 
enhanced for u>l — 2.6 eV, is similar to the one cor- 
responding to Si, S2 and S3 excitations as well as to 
the behavior of the new modes seen around 55, 59 and 
68 cm -1 . The results of perturbational analysis regard- 
ing energy scales and symmetries^, in particular of a Bi 
symmetric 5 = states at 45 cm -1 , argues for the mag- 
netic nature of the mode S3 , Refill. Just on the account 
of the similar resonance behavior of the S3 excitation 
with respect to the modes found above 50 cm -1 one can- 
not unambiguously identify the latter as magnetic bound 
states as well. The lack of energy shift in magnetic fields 
is consistent with such an interpretation, but the same 
would be true if they had a phononic origin. 

We discuss below the nature of the two light cou- 
pling mechanisms to magnetic excitations. For the set 
of modes discussed in Fig. 6c-d we propose that the cou- 
pling takes place via the spin-orbit interaction which can 
be written in an effective form as iJ™* oc (ej„ x e out )S t * t 
(Ref^). As it is found experimentally, in zero field 
this interaction Hamiltonian probes indeed excitations 
with A2 symmetry (the Tif, mode) and the calculated 
magnetic field dependent intensities for this and other 
modes in Fig. 6c-d is also in agreement with the ex- 
perimental results. The coupling to the Ti e and Tif 
modes from Fig. 3 can be possibly understood if we in- 
voke the usual effective spin interaction corresponding 
to the photon induced spin exchange process H™* oc 
J2<i,j>( e in r ij){eout*ij)Si ■ Sj. Here the sum runs over 
pairs of lattice sites, Si and Sj are the exchanged spins on 
sites i and j respectively, while is the vector connect- 
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ing these sites^. The explicit expression of this interac- 
tion for several polarizations in the 4-spin cluster approx- 
imation contains finite coupling in Bi and B2 channels 
for the triplet Ti e and Ti / states. This explains the pres- 
ence of the 25.6 cm -1 magnetic modes in all polarizations 
except (RR). 

The difference in the coupling strengths seen in Fig. 7 
is thus understandable because these two light coupling 
mechanisms given by 7JJ n * and H^ 1 * need not be simul- 
taneously in resonance with the same high energy ex- 
cited electronic states. The photon induced spin ex- 
change Hamiltonian, H™ 1 , is usually invoked in order 
to explain Raman active 5 = two-magnon type exci- 
tations in various magnetic systems 2 ^. SrCu2(BC>3)2 is 
an example where this Hamiltonian, in the presence of 
singlet-triplet mixing DM interactions, can be used to 
account for coupling to S = 1 states. We also note that 
in principle the photon induced spin exchange could pro- 
vide coupling to the 21.5 cm -1 5 = bound state be- 
low the spin gap. The reason this mode is not directly 
observed in our spectra for any of the two excitation en- 
ergies used in Fig. 6 is an open question. However, one 
possible explanation is that the Raman form factor for 
exciting pairs of k = magnons is vanishing as opposed 
to the case zone boundary modes. An example is the 
case of the Raman vertex calculated for the 2D square 
lattice within the spin-wave approximation 2 and using 
the Fleury-Loudon2£ interaction. Consequently, both the 
21.5 and the 28.9 cm -1 excitations could be attributed to 
5 = bound states but originating from different parts 
of the Brillouin zone and having substantially different 
binding energies. 

VII. SUMMARY 

We study by Raman scattering collective magnetic ex- 
citations in the spin-dimer compound SrCu2 (603)2 • Re- 



garding the one-triplet sector, we showed that by using 
a 4-spin cluster approximation and by including an addi- 
tional intra-dimer DM interaction we are able to explain 
the observed zero field symmetry selection rules and the 
rich behavior in magnetic fields. We are also able to ex- 
perimentally demonstrate the existence of two effective 
magnetic light scattering Hamiltonians which are respon- 
sible for their resonance behavior. The 4-spin approxima- 
tion fails to account for the excitations seen in the multi- 
particle magnetic sectors. In particular, the existence of 
a set of four modes below the onset of the two-triplet 
continuum (at 37.5, 40.8, 44.5 and 50.9 cm -1 in the A 2 
symmetry channel) shows that further theoretical analy- 
sis is required in order to understand the nature of these 
composite excitations. We suggest a possible explanation 
for the existence of a sub-gap collective mode in terms 
of a strongly bound singlet state which can be generated 
within the space of 4 nearest neighbor spin dimers. 

Experimental data in the energy range below 350 cm -1 
shows the existence of several quasi-degenerate phonons. 
General symmetry arguments suggest that these excita- 
tions, which are very close in energy, involve different 
vibrational patterns (in plane and c-axis atomic motions 
respectively). The failure of group theory to provide an 
understanding of this interesting behavior at a qualita- 
tive level calls for a detailed theoretical investigation of 
the vibrational modes. We also find a weak A 2 sym- 
metric phonon at 155.9 cm -1 . The energy and symme- 
try of this excitation suggest that this is a coupled spin- 
phonon excitation and can be evidence for the existence 
of magneto-elastic interactions in SrCu2(B03)2 • 
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